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, Abstract 

We investigate the construction of black holes and black strings in vacuum 
plane wave spacetimes using the method of matched asymptotic expansions. 
We find solutions of the linearised equations of motion in the asymptotic re- 
gion for a general source on a plane wave background. We observe that these 
solutions do not satisfy our previously defined conditions for being asymptot- 
ically plane wave. Hence, the space of asymptotically plane wave solutions is 
restricted. We consider the solution in the near region, treating the plane wave 
, ' as a perturbation of a black object, and find that there is a regular black string 

solution but no regular black hole solution. 
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Plane waves are of interest both from the point of view of classical gravity and in 
the context of considerations of holography in string theory. To a relativist, the 
plane waves are a rich class of exact solutions, which can be obtained as the result 
of applying the Penrose limit to an arbitrary spacetime. In string theory, they are 
of interest because the theory on the worldsheet admits a simple realisation, making 
^ explicit computations possible [I] . Secondly, certain maximally supersymmetric plane 

wave backgrounds [2] admit a dual field theory interpretation as a scaling limit of 
certain field theories [3j. This correspondence is obtained by taking a Penrose limit 
of the AdS/CFT correspondence. As a result, our understanding of holography in 
this case is rather indirect, and no holographic dictionary has yet been constructed 
for this duality. 

From both these points of view, the study of black hole solutions with plane wave 
asymptotics is clearly interesting. Black holes for other "simple" asymptotics, such 
as flat space or anti-de Sitter space, have long been known. The construction of 
black hole solutions with plane wave asymptotics would offer a new, rich family of 
black hole solutions, whose thermodynamics could exhibit interesting dependence on 
the asymptotic plane wave considered. For string theory, the black hole solutions in 
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the maximally supersymmetric plane wave backgrounds of [2] would presumably be 
related to finite-temperature excitations of the corresponding field theory, as in the 
AdS/CFT correspondence |3]. Consideration of such solutions could therefore cast 
interesting light on this poorly-understood duality 

Some exact solutions describing black strings in plane wave backgrounds have been 
obtained by applying solution generating transformations [3 EJ [7J, El E] . A review of 
this work and the structure of horizons and plane waves can be found in [TD] . However, 
such methods are available only in special cases, and a solution describing the simplest 
situation, a regular black hole or black string in a vacuum plane wave background, 
has not been obtained by these methods. Constructing solutions by directly solving 
the equations of motion is challenging. 

In this paper, we adopt the method of matched asymptotic expansions to find 
approximate solutions when the horizon size r + of the black hole or black string is 
small compared to the curvature scale [i~ x of the plane wave. This gives a separation 
of scales, which can be exploited to solve the equations of motion in the linearised 
approximation in separate regions, matching the solutions in an overlap region. Such 
methods have been successfully applied to the construction of caged black holes in 
Kaluza-Klein theory [TT], [T2"] and to construct black ring solutions in more than five 
spacetime dimensions [13] and in anti-de Sitter space [13] ■ These ideas have been 
further developed in [T5| ITS], where general extended black objects wrapping a sub- 
manifold in an arbitrary spacetime have been considered at leading order in the region 
far from the black object. 

We proceed in a similar way to these previous examples, first finding the metric 
far from the source (for r ^> r + ) by studying the linearised approximation to gravity 
with an appropriate delta-function source. The wave equation in the plane wave back- 
ground is rather complicated, so we focus on solving this problem in an intermediate 
region r + <C r /i -1 where the deviations from flat space due to both the source 
and the plane wave are small. 

Solving the equation in this regime, we find that simple dimensional analysis in- 
dicates that the solutions will violate the asymptotic boundary conditions proposed 
in [T7J as a definition of asymptotically plane wave spacetimes. In fact, the perturba- 
tion due to the delta-function source becomes large relative to the background metric 
at large distances. An explicit analysis in four and five dimensions shows that the 
terms violating these boundary conditions are indeed non-zero. These solutions thus 
appear not to be asymptotically plane wave; we will refer to them as black holes or 
black strings in plane wave backgrounds. The fact that the linearised solutions for a 
delta- function source violate the asymptotic boundary conditions suggests that as in 
AdS 2 [IB] and the Kerr/CFT correspondence [HI EHl EI] , the space of asymptotically 
plane wave spacetimes may be highly restricted. 

We then obtain the near horizon metric in the region r <C by solving the 
linearised Einstein equations on the background of the black object, treating the 
plane wave as a perturbation. For a black hole, we find that there is no linearised 
solution which is regular on the horizon. For the black string, we obtain a regular 
solution in the near region, and verify that it matches on to the solution in the 
intermediate region. 
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When solving the equations, we focus on vacuum plane waves in the lowest possible 
dimension for simplicity, but the method of matched asymptotic expansion is more 
general, and a similar analysis could be applied to construct black string solutions 
in any plane wave background of interest in arbitrary dimensions. We will remark 
on the extension to other waves and higher dimensions at appropriate points in the 
calculation. The calculation in the region r ^> r + is described in section [2], and the 
calculation in the region r <C /i -1 is described in section El We conclude with some 
remarks on the implications of our results in section HI 



2 Linearised solutions on a plane wave background 

We want to construct solutions corresponding to a black hole or black string of radius 
r + in a general vacuum plane wave background in D = d + 2 dimensions 

ds 2 = -dt 2 + dz 2 + fx aP (t + z)x a x p (dt + dz) 2 + 5 af3 dx a dx> 3 , (1) 

where x a , a = 1, . . . d are Cartesian coordinates on the transverse space. We will work 
in the parameter range r + <C yT , where we take the matrix n a p(t-\-z} characterising 
the wave to have a single characteristic scale \x for simplicity. The black object can 
then be treated as a small perturbation of the plane wave background for r ^> r + . In 
this region of the spacetime, the problem of constructing a black hole or black string 
solution thus reduces to solving the linearised Einstein's equations for a suitable source 
T^y. In transverse gauge, the linearised equations ar^j] 

□V = —lQnGT^y. (2) 

For a pointlike source, the relevant stress tensor is simply = MV^V v 5{x^ — x M (r)), 
where x M (r) is the particle's trajectory, = dx^/dr is the tangent to this trajectory, 
and M is the proper mass. For a black string solution, the stress tensor can be 
determined by linearising the vacuum black string solution in d + 2 dimensions, 



r d-2\ / r d-2 



ds 2 = -[l-^ dt 2 + dz 2 +ll--±^) dr' + r'dnU, (3) 



-i 



which gives the stress tensor in these coordinates as 



„d-2 „d-2 



The source is fixed to follow some appropriate trajectory in the plane wave back- 
ground. For a pointlike source, the appropriate trajectory is a timelike geodesic of the 
background spacetime. To obtain a stationary black hole solution, we should require 
this geodesic to be the orbit of a timelike Killing vector in the spacetime. This forces 
us to restrict to plane waves with a constant matrix \i a ${t + z) — fi a/ 3, so that the 



-'^Note that in our actual calculations, we will not assume the transverse traceless gauge, as it is 
more convenient to use the gauge freedom to fix particular components of the perturbation. 
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solution has a timelike Killing vector, and to consider the geodesic z = x a = 0, 
which is the unique geodesic trajectory which is also an orbit of the Killing vector. 
The appropriate source is then T tt = M5(z)S d (x a ), and the size of the black hole is 
r^T 1 oc M. 

For the black string, the equation of motion for a probe string is [22] 

V r = °> ( 5 ) 

where K p is the second fundamental tensor of the submanifold defining the em- 
bedding of the string worldvolume, and T^ v is the stress tensor of the source. We 
will consider embedding the black string along the submanifold x a = 0, which has 
K p = 0. As a result, there is no constraint on the form of the stress tensor. As 
for the black hole, we need to restrict to constant n a p(t + z) = n a p so that this sub- 
manifold is an orbit of the spacetime isometries, so that we can expect to obtain a 
stationary uniform black string solution. We can then use boosts in the t — z plane to 
choose the black string solution to be in its rest frame, setting T tz = 0, without loss 
of generality. The appropriate source is thus (Sj). We want to find a uniform black 
string solution, so the components of the stress tensor are assumed to be constants 
along the worldvolume. The blackfold equations of [16] are hence trivially satisfied. 

In each case, the problem thus reduces in principle to solving (J2]) on the plane 
wave background for an appropriate source. However, we do not have the Green's 
function for this differential equation in closed form, so we will content ourselves with 
studying this problem in the intermediate region r + <r< pT , where we can treat 
the plane wave itself as a small perturbation of flat space, and obtain the solution of 
t^j order by order in /i 2 r 2 . 



2.1 Dimensional analysis 

We first discuss the perturbation in general dimensions using a simple dimensional 
analysis argument. For the case of a point source, we find it convenient to rewrite 
the metric in spherical polar coordinates, introducing a radial coordinate 

r 2 = z 2 + 8 aP x a x^ (6) 

and defining coordinates 9 l on the S d at constant r. As in [23], we use a, b to denote 
coordinates on the two dimensional space spanned by r, t. By dimensional analysis, 
the form of the perturbation to first order in M and in /z 2 will be 

M (o) M/i 2 (1) 
M (o) M/i 2 (1) 

h- = JL h w + M£ h e>rgi) 

' k r D-5 U V ^ r D-7"'i3 \° » 



2 We can make this choice without loss of generality by translation invariance in z. 
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where h^} and hfy are dimensionless functions depending only on the angles 9 % . In 
fact, since the spherical symmetry is only broken by the plane wave, are constants, 
and the component on the sphere h\j will be proportional to the metric on the sphere 

7ij We will always work in a gauge where hf^ vanishes. Each addition of an i index 
raises the power of r by one because the coordinates on the sphere are written in 
terms of dimensionless angles. 

This simple dimensional analysis already indicates a significant issue: this per- 
turbation does not satisfy the boundary conditions introduced in [17]. There, it was 



assumed that components of the perturbation in the directions transverse to the wave 
would fall off at least as l/r D ~ 4 (corresponding to hy oc l/r D ~ 6 , because of the extra 
factors of r from writing the perturbation in polar coordinates), characteristic of a 
localised source in a flat spacetime. However, we find that the term resulting from the 
interaction with the wave must grow more quickly than this on dimensional grounds. 
When we think of the plane wave as a perturbation around flat space, the plane wave 
background introduces corrections which grow more quickly with r than the original 
leading-order response. 

Similarly, when we consider a black string source, it is convenient to write the 
metric in the directions transverse to the wave in polar coordinates, introducing a 
radial coordinate 

r 2 = 5 a(3 x a x 13 , (8) 

and introducing coordinates 9 l on the S^ 1 at constant r, z. In the string source case, 
a, b will denote coordinates in the three dimensional space spanned by t, r, z. Then 
to leading order in r + and /z 2 , the perturbation sourced by a black string will have 
the form 

h h - r+ h {0) i + ^ h (1) (e i ) 

nab ~ r D-A n ab + r D_ 6 a ab \P J> 

D— 4 D~ 4 2 

k ™ ~ ^E k ai + rD _ 7 h ai \P h K 9 ) 

D—A n-4,,2 

h V - ^6 h ij + r D-8 h i3 ^ h 

where h^} are constants and h^l are functions of the coordinates 9 l on the sphere 
only. Thus, as in the black hole case, the perturbation does not satisfy the boundary 
conditions introduced in [T7] . 

This is a significant issue because at least in low spacetime dimensions, the re- 
sulting perturbation actually grows more quickly with r than the background metric. 
In D = 4 for the black hole and D = 5 for the black string, the perturbation of the 
angular metric hij has a contribution that goes like r + /i 2 r 3 , which is growing faster 
than the background metric on the sphere, which goes like r 2 . Furthermore, what we 
have discussed so far is just the leading order correction in /z 2 . Higher order terms in 
fi 2 will come with additional powers of r. One might hope that when the problem is 
solved to all orders in /i 2 , the resulting behaviour could be under better control, but 
it is hard to see how such a cancellation between different orders could be arranged. 
We will see later in a particular example that this does not occur. 
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Thus, we are faced with the odd situation that the linearised field of a point source 
may become more important than the background, signalling a breakdown of the 
linearised approximation, far from the source itself. Thus, the solutions we construct 
should not be thought of as "asymptotically plane wave" black holes/strings, as the 
metric in the asymptotic regime is not close to the original plane wave metric. As 
a result, the analysis of [17] will not apply to these spacetimes, and in particular we 
do not expect that they will have finite action with respect to the action principle 
discussed there. 

One might hope that the terms which violate those boundary conditions which are 
allowed by dimensional analysis may actually vanish. This hope would be encouraged 
by the fact that the specific examples of plane wave black strings constructed in 
[51 El U\ satisfied the asymptotic boundary conditions of [T7] . However, the examples 
of [51 El El are special cases in that they are constructed by the Garfinkle-Vachaspati 
solution-generating transformation [21], and by construction can only differ from the 
seed solution in the metric components along the null direction. By contrast, the 
solution constructed in [9], which was obtained by a different method, has precisely 
the kinds of corrections that are predicted by this dimensional analysis argument. 

In the next two subsections, we will consider the solution of the linearised equa- 
tions of motion for the perturbation in detail for the lowest possible dimension for 
black hole and black string sources, and see in these particular examples that the 
terms which violate the asymptotic boundary conditions of [T7] do indeed appear. 
Thus, the approximate solutions we obtain for black holes and black strings in plane 
wave backgrounds are not asymptotically plane wave in the sense defined in |17j . 
Given the above dimensional analysis arguments and the results below, it seems rea- 
sonable to expect that this is the generic case, so that the space of asymptotically 
plane wave solutions is very limited. We will comment on this in the conclusions. 

2.2 Black hole 

Let us consider the perturbation sourced by a point source in the lowest possible 
dimension, D = 4, in detail. By a choice of coordinates, the most general four 
dimensional vacuum plane wave can be written as 

ds 2 wave = -dt 2 + dx 2 + dy 2 + dz 2 - [i 2 (x 2 - y 2 ) (dt + dz) 2 . (10) 

We rewrite this in spherical polars by defining 

z = rcos8, x = r sin 6* cos 0, y = r sin 6* sin 0, (11) 

so 

d s 2 wave = -dt 2 + dr 2 + r 2 (d9 2 + sin 2 6d(p 2 ) (12) 
— /i 2 r 2 sin 2 #(cos 2 (p — sin 2 (f))(dt + cosOdr — r sin 6d6) 2 . 

As in the above subsection, we can use dimensional analysis to fix the dependence 
of the perturbation on r. We can in fact determine the perturbation to zeroth order in 
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/i 2 by simply linearising the Schwarzschild solution, which gives h tt = h rr = — . This 
satisfies the linearised equations of motion for a delta-function point source, but not 
in the transverse traceless gauge which was assumed in writing (121) . In what follows, 
we will not assume the transverse traceless gauge, as it is more convenient to use the 
gauge freedom to fix some components of the perturbation. 

For the terms of first order in fi 2 , we can use the freedom to choose a gauge 
for the perturbation to set h^J and hfj to zero. Note that we have four gauge 
degrees of freedom but have only eliminated three components, hence we have one 
remaining degree of freedom which we will use later. We then make an ansatz for the 
dependence of the perturbation, and write our perturbation as 

hab = yh^ + M^r(cos 2 <P-sm 2 <P)h^(e), 

h a6 = MA 2 (cos 2 0-sin 2 0)/i2/(0), (13) 
= M y V 3 '(cos 2 <f>- sin 2 <j>)h§ '(0), 

where the non-zero components of are f% = 2, h$ = 2, and the non-zero 
components of h$(9) are h$(9), h ( t 1 r ) (9), h$(9), h$(9), h^(9), h$(9) and h${6). 

We now want to substitute this ansatz into the linearised Einstein equations and 
solve for the undetermined functions h$(9), requiring regularity on the sphere. In 
an arbitrary gauge, the linearised Einstein equations for r ^ are 

4- = 7^ p<T (v P v M /w + v p v,/v - v M v„v - v p v ff v) = °- ( 14 ) 

Substituting our ansatz, these equations become (where primes denote derivatives 
with respect to 9) 

- sin 2 9hS ) "(9) - sin 9 cos 9h$'{9) + 2(cos 2 + l)h { ^(9) (15) 
-6 cos 6 9-2 cos 4 9 + 22 cos 2 9 - 14 = 0, 

- sin 2 9h < £ ) " (9) -sin9cos9h ( t l ) \9) + 2 sin 2 0/^ } '(0) +4h${6) (16) 
+2 sin 9 cos 6h$ (9) - 8 cos 5 + 16 cos 3 9 - 8 cos 9 = 0, 

- sin 2 9h$" {9) -sin 9 cos 9h$' (9) + 4 sin 2 Ohty' (0) + 2 (3 - cos 2 0) hg> (0) (17) 
-2 sin 2 0h$ (0) - 2 sin 2 Ohfy (0) + 10 cos 6 - 26 cos 4 + 22 cos 2 0-6 = 0, 

sin 2 9h^'{9) - sin 2 9h { ${9) - sin cos 9h$ (0) + sin cos 9hfJ (0) (18) 
+2(cos 2 + 1)^(0) -4sin0(cos 5 0-2cos 3 + cos0) = 0, 

sin 2 0h$\6) + 2(cos 2 + l)h% ] (0) + 2 sin 0(cos 4 0-4 cos 2 9 + 3) = 0, (19) 
/$'(0) - cot 91$ (0) - h§ (0) + /i£> (0) - cos 4 + 2 cos 2 - 1 = 0, (20) 
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hf/{9) -cot 0/1$ (0) + h\l\e) + 2 cos 0(1 -cos 2 0) = 0, (21) 

-sin 9(h$'(9) -h$'{9)) +cos 0(h$ (0) -h$ (0)) +2 sin 0/4 1 , ) (0)+2cos0(cos 4 0-l) = 0, 

(22) 

- sin 2 6hf;(6) + sin 2 9h$"(9) + sin cos 9{h§'{9) - 2h${9)) 

- sin 2 0/i$"(0)6 sin 2 0/i$'(0) + -5 sin 2 9h$(9) + 3 sin 2 6h$(6) (23) 
- sin 2 0/i$ (9) + sin 2 (0) + 2 cos 2 0(cos 4 5 + 3 cos 2 9 - 5) + 2 = 0, 

sin 2 <i>) + sin0 cos 0(/i«'(0) - h$'(9) + 2h${0) - h$>'(0)) 
+ cos 2 0(3/$ (9) - 3h$ (9) + h$ (0)) - 2 sin 2 9h$'{9) + 3 sin 2 0h% (9) (24) 
-7hg) {9) + 3/4 1J (9) + 2 cos 2 0(3 cos 4 9-7 cos 2 + 9) - 10 = 0. 

We have a system of ten equations in seven unknown functions (in fact, there will 
be only six unknown functions once we have made use of the one remaining degree 
of gauge freedom) and so it seems our system is over-constrained. We find, however, 
that there are only six independant equations and hence that our system is in fact 
well defined. It is convenient to subtract a multiple of ( l20l) from ( fi~7l) to simplify it to 

- sin 2 9hg}" (0) - sin^ cos 6h®'(6) + 2(1 + cos 2 6)h$(0) + 8 sin 9 cos 6h$ (9) (25) 
+2 sin 2 9h$ (9) - 2 sin 2 (0) + 10 cos 2 0-14 cos 4 + 6 cos 6 0-2 = 0. 

By using combinations of f[2"Uj) . (j2"2l . (1251) and their derivatives it is possible to reduce 
(fT5"j) to an algebraic equation 

2 sin cos 0/i$ (0) + 2 sin 2 0/i$ (0) + 3/4? (0) - 5/i« (0) (26) 
+2 cos 2 9h$ (0) + 2(- cos 6 + 4 cos 4 + cos 2 - 4) = 0. 

We find we can write (fB]l. (jT5]). (|25j) and (Elj) as linear combinations of (TT9"j) . (l2Uj) . 
f[2"Tj) . f[2"2"j) . (125]) and f[2"6l) and hence that these equations are not independant. We 
now see that a convenient choice of gauge is one in which h$(6) = 0. In this eauee, 
the solution which is regular on the sphere is 

h^dx^dx" = ^-dt 2 + ^-dr 2 + Mfi 2 r sin 2 0(cos 2 <p ~ sin 2 4>) x 

(4 - i sin 2 0)dt 2 + 4 cos 0dtdr - 4r sin 0rftrf0 (27) 

+^ sin 2 9dr 2 - jjr 2 sin 2 0(d0 2 + sin 2 9d<f) 2 ) . 

We note that as stated earlier, the regular solution for the terms of first order in 
fi 2 has non-zero components on the sphere which grow faster than the background 
metric on the sphere. These solutions are hence not asymptotically plane wave. While 
this leading order term would not grow faster than the background metric in higher 
dimensions, higher order terms in /i 2 will in principle do so. 



8 



2.3 Black string 



We now consider the perturbation for a black string source in the lowest possible 
dimension, which is D = 5 for the black string. The most general vacuum plane wave 
solution in five dimensions is 

ds 2 wave = -dt 2 +dx 2 +dy 2 +dz 2 +dw 2 -/j, 2 (a(x 2 +y 2 -2w 2 )+(3(x 2 -y 2 ))(dt+dz) 2 ; (28) 

note that there is a two-parameter family of plane wave solutions here. We rewrite 
this in spherical polars in the directions transverse to the wave by writing 

x = r sin 9 cos 0, y = r sin 9 sin 0, w = rcos9, (29) 

so 

ds 2 wave = -dt 2 + dz 2 + dr 2 + r 2 (d9 2 + sin 2 9d<p 2 ) (30) 
-H 2 {ar 2 (l - 3 cos 2 9) + (3r 2 sin 2 #(cos 2 - sin 2 <f>))(dt + dzf. 

As in the previous subsection, we can determine the perturbation to zeroth order 
in fi 2 by simply linearising the Schwarzschild black string solution ([3]), which gives 
h tt = h rr = 2 ^-. We will again find it convenient to fix the gauge by choosing some 
components of the perturbation to vanish at each order in /i 2 . For the terms of first 
order in fi 2 , we note that the background has an invariance under t — > —t, z — > —z 
which is not broken by the source, so the h tfl , h Zfl components for fi ^ t,z will 
automatically vanish. 

At first order in /i 2 , we can treat the two different components of the plane wave 
separately. We first consider the first-order terms in the perturbation associated to 
a. Let us therefore set a = 1 and (3 = in the plane wave background ( 1281) . There is 
then a translation invariance in and a symmetry under — > —0, which imply that 
hfa vanish for \i ^ 0. We will make a choice of gauge to set hrJ and /i^ to zero. 
This gauge choice proves to be convenient for comparing to the solution in the near 
region to be obtained later. The form of the perturbation is then 

h ae = Mfi 2 r 2 h^(9), (31) 
fHj = Mfi 2 r 3 h^(9), 

where the non-zero components of h^ are h# = 2, h$ = 2, and the non-zero 
components of h$(9) are h$(6), $(6), h^{9), h$(6) and h%(6). 

We now want to substitute this ansatz into the linearised Einstein equations and 
solve for the undetermined functions h^}{9), requiring regularity on the sphere. In 
an arbitrary gauge, the linearised Einstein equations for r ^ are 

R$ = ^(VpVA, + v p v,/v - v M v, v - v p v ff v) = °- ( 32 ) 
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-3cos 2 0) = 0, 


(33) 


-3cos 2 0) = 0, 


(34) 


_3cos 2 0) = 0, 


(35) 


) = o, 


(36) 


sin 2 = 0, 


(37) 




(38) 



Substituting our ansatz, these equations become 

d 2 h$(9) + cot 9d e h$ (9) + 2h$(9) + 1 

dgh^Jie) + cot 0^/4? (0) + 2h S( d ) + 12(1 - 3 cos 
%f>f$ (0) + cot Bd e h®{6) + 2h®(6) + 8( 
/iS(0)+/ i 2(0) + 2(l-3cos 2 

tan Odeh^l (0) - /»« (0) + fcg (0) + 6 sin 2 = 

9 2 / i ;; ) (0) - %h%{0) - %hs>(0) + cot 0( 

- - M<8(0) - h%{e) + 12sin 2 - 2(1 - 3cos 2 0) = 0, 

dlh^{9) + cot B{deh®(9) + d e h^(9) - - rf^)) (39) 

+3/4J(0) + 3/^(0) - ^(0) + 45(0) + 2(1-3 cos 2 0) = 0. 

We first solve equations (1331) . ( 1341 and f l35l) for /i# (0), h®(6) and h®{9) respectively. 
We then solve for h^(d) and h^l(9) using equations (TTTT) and (1371) . It is easy to 
verify that these solutions satisfy ( 1381) and (1391) . Keeping only the regular part of the 
solution, we find 

h$(9) = 4(1 -3 cos 2 0), /£ } (0) = 3(1 -3cos 2 0), h®{6) = 2(1 - 3 cos 2 0), (40) 

hfj{6) = -(1 -3cos 2 0), /i$(0) = -sin 2 0(1 -3cos 2 0). (41) 

As in the black hole case, we see that terms that grow faster than the background 
metric at large r do indeed occur. 

It turns out that for this background, the linearised equations of motion can be 
solved exactly by including one further term at next order in fi 2 . If we take 

h ab = yh^ + M^rh ( 2(9) + M^r 3 h < S(9), 

h ae = M^ 2 r 2 h^(0), (42) 
h i3 = Mfi 2 r 3 hff(9), 

with hjO!v and as given above, and 

>4 a) = 4 2) = h$ = 1(3- 30 cos 2 + 27 cos 4 0), (43) 

this will solve the equations to linear order in M but to all orders in /i 2 . This gives 
an approximation valid in the full far region r ^> M, demonstrating that the bad 
asymptotic behaviour of this solution is not resolved at higher order in /x 2 . 
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We now consider briefly the similar analysis for the other independent component, 
setting a = and (3 = 1 in the plane wave background (1251) . The <f> dependence in 
this background restricts our ability to simplify the form of the solution by general 
arguments, but the results from the previous case suggest we take an ansatz of the 
form 

h ab = —h$ + M/i 2 r sin 2 9 (cos 2 0- sin 2 </>)h$ , 
h aB = M[i 2 r 2 sin 2 9 (cos 2 0- sin 2 <j>)b$ , 

h ee = M[i 2 r 3 sin 2 9{cos 2 <f>- sin 2 <f))h$, (44) 
h H = M/j, 2 r 3 sin 4 #(cos 2 0- sin 2 <f>)h^, 

assuming the angular dependence at first order in fi 2 will reproduce the angular 
dependence of the background plane wave. The non-zero components of are 
hft = 2, h^J = 2, and we assume the h$ above are constants. We find that we can 
solve the linearised equations of motion to first order in /i 2 for this ansatz by setting 
_ 4 _ 3 iS}) - o ?,(i) - _i h {1) - -1 

We can summarise these results in a more invariant fashion by saying that for a 
plane wave background of the form 

ds 2 wave = -dt 2 + dz 2 + dr 2 + r 2 {dd 2 + sin 2 6d<p 2 ) - fi 2 f{6, <j>){dt + dz) 2 , (45) 

a solution of the linearized equations of motion for a black string source, to linear 
order in /i 2 , is 

htt = — + AMfi 2 rf(6,<p), (46) 
r 

h tz = 3Mfi 2 rf(9,^, 

h zz = 2M/i 2 r/(M), 
2M 



h 



rr 



r 

,2^3 . 



hee = -M^r'fiO,^), 
h H = -M^ 2 r 3 sin 2 9f{9,(j)). 

We would expect that this generalises straightforwardly to higher dimensions. As 
in the black hole case, this demonstrates that these solutions are not asymptotically 
plane wave, as the perturbation is large compared to the background metric far from 
the source. 



3 Near region analysis 

Having explored the behaviour in the intermediate region, where we can use a lin- 
earised approximation about the plane wave background, we now turn to the analysis 
in the region r ^ fi^ 1 near the black hole or black string. In this region we can 
treat the plane wave as a small perturbation of the black object, and the problem 
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reduces to linearised perturbations on the black hole or black string background, with 
boundary conditions at large distances determined from the previous solution in the 
intermediate region and a boundary condition at the horizon determined by requir- 
ing regularity of the perturbed solution there. We will find that there is no regular 
solution in the black hole case. For the black string, we find a regular solution which 
matches on to the solution we discussed above in the intermediate region. We will 
focus on the analysis for the black hole in four dimensions and the black string in five 
dimensions, as in the previous section, but the same techniques can easily be applied 
in higher dimensions. We will comment briefly on the extension of the analysis to 
higher dimensions for the black hole case. 

3.1 Black hole 

We first study the near region of the black hole, treating the plane wave as a pertur- 
bation. We will do the analysis in the lowest possible dimension, D = 4, even though 
there is a simple symmetry argument that no regular solution exists in this case. The 
calculation is simplest in this dimension, and it serves to illustrate the method of 
calculation, which will be very similar in higher dimensions. 

Take the Schwarzschild black hole solution in four dimensions, 

dv 2 

ds 2 = -f(r)dt 2 + -7TT + r 2 (d9 2 + sin 2 8d<j) 2 ), (47) 
fi r ) 

with f(r) = l — 2M/r. We want to find a solution of the source- free linearised vacuum 
equations on this background which asymptotically approaches the four-dimensional 
plane wave ffTUj) . This implies that we want a perturbation with asymptotic 
boundary conditions 

2 2 

lim huvdx»dx v = -— sin 2 6(e 2i * - e - 2i ^)(dt + cos 6dr - r sin OdO) 2 + ..., (48) 

r^oo 2 

where the . . . denotes terms going like [i 2 M n for n^O. These terms are suppressed 
relative to the leading term because dimensional analysis tells us the mass will always 
appear in the combination M/r. At linear order in M the subleading terms at large 
r should match onto the results of the analysis in the intermediate region obtained 
in the previous section. 

At the horizon, the boundary condition is that the solution be regular there. Since 
the background metric is not regular at the horizon in the Schwarzschild coordinate 
system we are using, this condition is most easily applied by writing the perturbation 
in an orthonormal frame. A suitable frame is = a/ f(r)dt, e*- 1 -* = f{r)~ l / 2 dr, 
e^ 2 ' = rd6, = r sin Od<p. Requiring that the components of the perturbation in 
the orthonormal frame are regular at the horizon implies that we must require that 
as r -> 2M, 

h u ~ (r - 2M) , V ~ (r - 2M) l ' 2 for ^ t, (49) 
hrr ~ (r-2M)"\ hr? ~ (r - 2My 1/2 for /i ^ r, (50) 
hij ~ (r - 2M)°. (51) 
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These conditions can also be derived by requiring finiteness of h^ v in a coordinate 
system which is well-behaved at r = 2M, such as Kruskal coordinates. 

Matching the leading term written in (|48p and imposing regularity at the horizon 
should determine the solution of the perturbation equations uniquely. In fact, as we 
mentioned above, we will find that there is no solution of the linearised perturbation 
equations that satisfies these two boundary conditions. 

For the black hole case, the analysis of the components on the sphere is sufficiently 
complicated that it is useful to exploit the results of [23J on the spherical harmonic 
decomposition for perturbations of Schwarzschild and rewrite the linearised equations 
of motion in terms of gauge-invariant variables with respect to coordinate transfor- 
mations on the sphere. We therefore want to convert (j4"81 into boundary conditions 
for their gauge-invariant perturbations. Let a,b — t,r and i,j = 6,<p. Then we have 
boundary conditions which are scalars h a b, vectors h a i, and a tensor hij, for which 
the boundary condition only has an hgg component. Following [23J we expand the 
perturbation in terms of harmonics on S 2 : the scalar harmonics 

03 = -1(1 + 1)3, 1 = 0,1,2,..., (52) 

the vector harmonics 

nVi=(-l{l + l) + l)Vi, 1 = 1,2,3,..., (53) 

with DiV 1 = 0, and the transverse traceless tensor harmonics 

nr y =(-Z(/ + l)+2)T fi , 1 = 2,3,4,..., (54) 

with DiTj = 0, T- = 0. There are, however, no pure tensor harmonics T^- on S 2 . We 
use the notation □ = DiD 1 for the d'Alembertian operator on S 2 , where Di is the 
covariant derivative with respect to the metric 7^- on the unit two-sphere. 
In terms of these harmonics, the scalar components of the perturbation are 

h ab = Y,fabST- (55) 
1. m 

Note that here and hereafter we will omit the I, m indices on the coefficients f a b or 
equivalent in the general relations like this for brevity. The vector perturbations are 
decomposed into their scalar derived and pure vector components h ai = + h^, 
where 

*4 = rX;/a(~ASr), (56) 

l,m 

where k 2 = 1(1 + 1), and 

hl = rJ2fa(Vrh (57) 

l,m 

Similarly, the tensor part of the perturbation is decomposed into scalar derived, vector 
derived and pure tensor components = hfj + hfj + hfj, where 

fcg = 2r 2 YjLHirtijSr + HtSij), (58) 

l,m 
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where S t , = ^D^SJ" + f Sf 1 , 

hl = 2r*Y,H V L V lv (59) 

l,m 

where V^- = -^(AVj + DjVi) with fcy = /(/ + 1) - 1, and 

hl = 2r^H T L T tr (60) 

There are, however, no pure tensor harmonics on S 2 . 

Thus, to determine the boundary conditions for the gauge invariant variables, we 
must apply this expansion to f HHl) and find the asymptotic values for the unknown 
expansion coefficients. For scalar perturbations this is straightforward. Substituting 
([48]) into (155]) we are able to read off that 

lim (fu)f = lim (f rr )f = (61) 

r^oo Z r- -*oc 14 

lim (ftr)¥ = - A 2 , lim (f rr )f = (62) 

r— >oo r— >oo 14 

We now turn our attention to the vector perturbations. Since D l Vi = we have 
= 0, so 

DX^D'hi^rJ^faSr, (63) 

l,m 

where we have used D % DiS = —k 2 S. Explicit computation gives us the boundary 
conditions for the scalar derived vector coefficients, 

lim(/ r )f = -^, lim (/ t ) 3 ±2 = 2/A 2 , hm(/ r )f = ^. (64) 

r— >oo i r— *oo r— >oo 14 

To find the pure vector coefficients we write 

hi = h ai - h s ai = h at + rJ2 fa^S? 1 = r £ f a ffi™),. (65) 

l,m l,m 

Again, by explicit computation we find, 

Iim(/ t ^ = ^, lim (/ ;) 3 ±2 = ^. (66) 
Finally we consider the tensor perturbations. We can write 

/i\ = (^)Wr 2 ^iW, (67) 

l,m 

where we have used D % Vi = 0, T\ = 0, S\ = and Yi — 2- This allows us to easily 
show that 

^(H L )f = - 3 -^f, hm(^)f = ^. (68) 
14 



To find the scalar derived transverse modes we will need the following results, 

D'&Vij = 0, (69) 

D i D i Sij = ^^-S, (70) 

which are proved in an appendix. Using the above results along with D l Tij = 0, we 
find 

D'LPhij = D'Dfhfj (71) 

= 2r 2S ^{-k 2 H L S] n + H T D i D j S ij ) 
i, m 

= 2r 2 J2(-k 2 H L + H T i -^ Y ^)S. 

l,m 

We can now show that 

lim (if T )f = \un(H T )f = - 5 J^. (72) 

r— >oo i r—*oo \_Z • I 

To find the vector derived transverse modes we will use the identities 

D% = -^# 2 - 2)DjS, (73) 



and 



D% = -^r(#-W, (74) 



which we also prove in an appendix. Since D l Tij = 0, we have 



and using the results above we can write this as 



D%j = D% + D%, (75) 



D { hij = 2r 2 ^(H L ~^(k 2 - 2)H T )D 3 S + 2r 2 £ H^(k 2 v - l)V 3 . (76) 

l,m l,m V 

We are now able to show that 

hm(i^ 3 ±2 = -^. (77) 

Using Maple we find that = hfj + hy, so there are no pure tensor perturbations 
as expected. 

We now want to translate this into boundary conditions for the gauge-invariant 
variables introduced in [23J. For vector perturbations the gauge-invariant variable is 

F < = fa + W D ^T- (78) 
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For I = 2, lim r ^ 00 (/ 4 l/ )f = so lim^ F t = ^f, and we have F a = r - l t ab D\r<5>) 
so the boundary condition for the vector master function is 



lim C2 = (79) 



For / = 3, lim™(/ r v )f = ^ and linv^tf^ 2 = so F Q = 0: this mode 

is pure gauge. This is as we might expect: the r 2 behaviour of the plane wave is 
typical of an I = 2 spherical harmonic, so the higher I modes that seem to appear 
in our decomposition of the mode in terms of spherical harmonics ought to be pure 
gauge. 

For scalar perturbations, the gauge-invariant variables are [23] 

F = H L + l -H T + -{D a r)X a (80) 
2 r 

F ab = f ab + D a X b + D b X a (81) 

with 



The master variable $ is 



with 



X a = p(/„ + rD a H T ). (82) 



^2Z- r (X + >0 (83) 



X = F*-2F (84) 
Y = F;-2F (85) 
Z = 0. 



For / = 2 perturbations direct substitution gives us lim^ocX = /i r , Y = 0, Z = 
hence the boundary condition on $ is 



lim $f =2 = -^t (87) 

For the / = 3 and / = 4 modes we find the gauge-invariant variables F and F ab 
are zero, so these modes are pure gauge as expected. Thus, we are left with two 
non-trivial modes, the / = 2 scalar and the I = 2 vector modes. 

Having established which modes are non-zero and their boundary conditions we 
consider the bulk solution. For the vector mode the equation for the master field is 
[23] 

d r ((l - — )d r $) - -Ml + 1) - 3 • — ]$ = 0. (88) 

iy /yZ 'X' 

The boundary condition is lim^oo $^1 2 = ^5-, therefore we set $ = r 3 ip. This allows 
us to reduce the master equation (IBB]) to 

d r {r\l - —)d r iP) = 0. (89) 
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which has solution 



\ 1 1 1 1 , . 2M.\ , , , 

4> = a[ - + — + — + — + rln(l ))+b. (90) 

y V 8Mr 12MV 3 16M 3 r 2 16M 4 r 32M 5 v r ' ) K ' 

Solutions with o ^ are clearly not regular at r = 2M, therefore the solution for 
the vector master field is $ y = br 3 . The boundary condition at large r then requires 
b = However, the boundary condition at the horizon (H91 requires that hu and hu 
vanish at the horizon. This implies that f{ and hence FY also vanish at the horizon. 
Finally F l = r~ l D r (r&) implies that $ too must vanish at the horizon, which would 
require 6 = 0. Hence, there is no solution which satisfies the boundary conditions at 
both the horizon and infinity. 

Thus, there is no regular solution describing a four- dimensional black hole in the 
plane wave background ( flOl) . In fact, this is not a surprising result in four dimensions; 
the rigidity theorem [25] shows that regular black holes must be static or stationary 
axisymmetric, and the plane wave ( flOl) is not static and does not preserve a U(l) 
symmetry. Thus, the plane wave perturbation breaks too many of the symmetries of 
the black hole for a regular deformed black hole solution to be possible. 

One might hope to avoid this problem by considering a non-vacuum plane wave 
solution. We can for example consider in four dimensions the electromagnetic plane 
wave 

ds 2 wave = -dt 2 + dx 2 + dy 2 + dz 2 - [i 2 (x 2 + y 2 )(dt + dz) 2 (91) 
supported by the electric flux 

F = 2n(dt + dz) A dx. (92) 

This is also interesting as a simplified model of the maximally supersymmetric plane 
wave of [26J. Here, the metric perturbation preserves a U(l) symmetry, but this is 
broken by the gauge field, and as a result, we again do not expect to find a regular 
black hole solution. In this case, the problem is that the equation of motion for 
the gauge field on the Schwarzschild black hole background has no solution which is 
regular on the horizon and satisfies the boundary condition at large r. 

If we consider the situation in higher dimensions, the above rigidity argument does 
not apply, but there is still no regular solution. Take for example a six-dimensional 
Schwarzschild black hole and add as a perturbation the six-dimensional vacuum plane 
wave 

ds 2 wave = -dt 2 + dv 2 + dw 2 + dx 2 + dy 2 + dz 2 - /i 2 (f 2 + w 2 - x 2 - y 2 )(dt + dz) 2 . (93) 

This clearly preserves two U(l) isometries, in the x — y and v — w planes. However, 
if we rewrite this in spherical polars, there is again an I = 2 vector part to the 
perturbation in the decomposition into spherical harmonics. The analysis is very 
similar to the above four- dimensional case, and it is not possible to find a solution for 
the vector part of the perturbation that satisfies the plane wave boundary conditions 
at large distances and the regularity condition on the event horizon. In this case, the 
plane wave preserves two U(l) isometries on the S 14 surrounding the black hole, so 
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the above argument does not apply; a regular deformed black hole solution would not 
violate the conditions of [21] • This problem seems to be very general. In all cases 
we have explored in the vacuum Einstein equations, the plane wave has a vector 
part in the spherical harmonic decomposition, and it is not possible to find a regular 
perturbation of the black hole which satisfies the plane wave boundary condition. It 
would be interesting to understand the physical origins of this restriction further. 

3.2 Black string 

We next study the near horizon region of the black string, treating the plane wave as 
a perturbation. The background is the five-dimensional black string solution 



with f(r) = l—2M/r. We want to find a solution of the source- free linearised vacuum 
equations on this background which asymptotically approaches the five- dimensional 
plane wave ([28]) . This implies that we want a perturbation h^ v with asymptotic 
boundary conditions 



lim h^ u dx^dx v = -/iV[a(l-3 cos 2 9)+f3 sin 2 #(cos 2 0-sin 2 0)}(dt+dz) 2 + . . . , (95) 



where the . . . denotes terms going like fi 2 M n for These terms are suppressed 

relative to the leading term because dimensional analysis tells us the mass will always 
appear in the combination M/r. 

As in the analysis in the intermediate region, we will deal with the a and /3 
components separately. It will turn out that the analysis is identical in these two 
cases. In terms of the spherical harmonic analysis on the two-sphere, these are scalar- 
type perturbations, which excite the I = 2,m = and Z = 2,m = 2 harmonic modes 
respectively. In the linearised theory, we can assume that the perturbation has only 
these modes turned on. Since only scalar-type modes are excited, the analysis on the 
sphere is fairly simple, and we will follow he similar analysis by Emparan et al [13J. 

The boundary conditions, and hence the perturbation, are invariant under simul- 
taneously taking t — > — t, z — > — z, so the only modes we need to consider are h tt , 
htzi h zz , h rr , and the longitudinal and transverse scalar-derived perturbations on the 
sphere. 

We first consider only the / = 2,m = perturbation (we set /3 = 0). Assuming 
that only this spherical harmonic is excited, we can write the perturbation as 



htt = oc(l — 3 cos 2 9)a(r), h tz = a(l —3 cos 2 9)b(r), h zz = a(l — 3 cos 2 0)c(r), (96) 




(94) 



r 



oo 




(97) 




(98) 
(99) 



h H = ar 2 sin 2 #[(1-3 cos 2 0)g(r) + 3 sin 2 6h(r)]. 
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Note that g(r) is the coefficient of the longitudinal mode on the sphere, and h(r) is 
the coefficient of the transverse mode on the sphere. As in [TB] , there is a remaining 
coordinate freedom, under 

r + 7(r)(l -3cos 2 fl), 9 -> + 6/3(r) cos ^sin^, (100) 

with 

W = ~ , 7(7 9 ) / .. v 7(2M) = 0. (101) 
r(r — 2M) 

Similarly, for the I — 2,m — 2 pertubation (obtained by setting a = 0), we define 
ha = 13 sin 2 fl(cos 2 0- sin 2 <f>)a(r), h tz = (3 sin 2 fl(cos 2 - sin 2 0)6(r), (102) 

/i« = /3 sin 2 #(cos 2 - sin 2 0)c(r), (103) 

sin 2 fl(cos 2 0-sin 2 0) 
rr P (l-2M/r) I[h [ } 

he<p = fir 2 sin 9 cos # sin cos 4>h(r), (105) 

= /5r 2 [sin 2 fl(cos 2 - sin 2 $)g{r) - (cos 2 + l)(cos 2 - sin 2 <j>)h(r)], (106) 

^ = /5r 2 sin 2 #[sin 2 fl(cos 2 - sin 2 <f>)g(r) + (cos 2 + l)(cos 2 - sin 2 </>)h(r)]. (107) 

Now we have remaining coordinate freedom under 

r — > r + 7(r) sin 2 (9(cos 2 — sin 2 0), (108) 

-> # + 2/3(r)sin#cos#(cos 2 0-sin 2 0), (109) 
— > — 4/3 (r) sin 2 0cos0sin0, (HO) 

with 

We find both coordinate transformations produce identical shifts 

2M / , 2M j(r) \ 

7(r), / (r) - / (r) + ^2 T ' - — -^j , (112) 

2 

#(r) -> #(r) + -7(r) - 6/3(r), h(r) -»• h(r) + 2/3(r), (113) 
r 

while 6(r) and c(r) are unchanged. 

We want to consider combinations which are invariant under these coordinate 
transformations. B = b(r) and C = c(r) are already invariant. We define in addition 

A = a(r) + —(g(r) + 3h(r)), (114) 
r 

P=f{T) -£ lrWT)+3Hr)))+ Mm^m , (11B ) 
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dh g(r) + 3ft(r) 

dr (r-2M) v ; 



Note that in this section, primes denote derivatives with respect to r. As in [13J, 
the constant part of h(r) can be fixed using the constant part of /3{r). Using the 
gauge-invariant combinations basically amounts to setting g(r) = —3h(r), which 
can be achieved for r 7^ 2M by an appropriate choice of gauge. Because of the 
boundary condition in fillip . g(2M) + 3h(2M) is gauge- invariant. It will however 
not be determined by solving the equations of motion for the above gauge-invariant 
variables, and will have to be separately specified. It will turn out to be determined 
by requiring regularity of the solution at the horizon. 

For either a = or (3 = 0, substituting into the linearised Einstein equations 
gives the same system of equations for the unknown functions A, B, C, F, H' (keeping 
terms up to 0(/z 2 )), 

oc r 2 (r -2M) 2 A" + r(r -2M)(2r -5M)A' - M(r -2M) 2 C (117) 
-(6r(r - 2M) - 2M 2 )A + M(r - 2M) 2 F' + 6M(r - 2M) 2 H', 

B!£ oc r(r - 2M)B" + 2(r - 2M)B' - QB, (118) 
#2 oc r(r - 2M)C" + 2(r - M)C - QC, (119) 

Rg) oc r 2 (r - 2M) 2 A" - rM(r - 2M)A' + 2M(2r - 3M)A (120) 
-r(r - 2M fC" - M(r - 2M) 2 C + (2r - 3M)(r - 2M) 2 F' 
+6(r - 2M) 2 F + 6r(r - 2M) 3 H" + 6(2r - 3M)(r - 2M) 2 H', 

oc -r 2 (r-2M)A' + r(r-M)v4 + r(r-2M) 2 C"- (r - 2M) 2 C (121) 
-(r - 2M)(r - M)F - r(r - 2M) 2 H', 



R ee + -T27{ R u K r ( r ~ 2M ) A ' ~ ( 3r + 2M ) A ~ ( r ~ 2M) 2 C + 3(r - 2M)C 



sin^ 



+ (r - 2M) 2 F' + 5(r - 2M)F + 3r(r - 2M) 2 H" (122) 
+6(2r - 3M)(r - 2M)F', 



Res 1 -T^ r{ 11 oc -rA+(r-2M)C+(r-2M)F+r(r-2M) 2 i^ // +2(r-M)(r-2M)iJ , . 

sin 6* w 

(123) 

In fact it is easy to show that the linearised Einstein equations must be the same for 
both modes. The perturbation involves some I = 2 scalar harmonic, let's call it S , 
so 

h ab = f ab (r)S, h ai = f a (r)V t S, ha = f{r)Sg l3 + f(r)V i V i 5', (124) 

where i,j are coordinates on the two-sphere and a,b = t,r,z. Then the first order 
Ricci tensor constructed from the second covariant derivatives of will also depend 
on angular coordinates only through S and its derivatives. Using ViV l S = —65 and 
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the fact that the sphere is an Einstein space, so Rij = gij, one can eliminate extra 
derivatives of S, to leave us with 

R ( 2 = e ab (r)S } R$ = e a (r)V^, R$ = e(r)S 9ij + e'(r)V,V^. (125) 

Hence, the resulting equations e a b(r) = e tt (r) = e(r) = e'(r) = are independant of 
whether S is in the m = or m = 2 mode. Thus, solving the equations (1117111231) 
will give us the general solution for the perturbation in the near-horizon region for 
both modes. 

The boundary conditions at large r imply that at order M°, a(r),b(r),c(r) — > 
—/j, 2 r 2 , and f(r), g(r),h(r) have no /i 2 M° term. This implies that 

A, B, C -/i 2 r 2 , (126) 

and F and H' have no /j, 2 M° term. Regularity at the horizon requires a(r) oc (r — 2M), 
b(r) oc \/r — 2M, and the other functions c(r), f(r),g(r) and /i(r) are required to be 
finite there. In terms of the gauge-invariant combinations, these boundary conditions 
are best expressed in terms of the alternative combinations 

A = A (r — 2M)H', F = F — MR' . (127) 

r 

The conditions for regularity at the horizon are then that A — > 0, F is finite, and H' 
is allowed to diverge like (r — 2M) _1 . 

We now want to solve this system of equations. We see that there are two de- 
coupled equations, (11181) and fll 19j) . The solutions of these satisfying our boundary 
conditions are 

B(r) = -/i 2 (r - M)(r - 2M) (128) 

and 

C(r) = -fi 2 (r 2 -2Mr + -M 2 ). (129) 

3 

It is also convenient to subtract a multiple of (11191) from fll 20[) to simplify it to 

= r 2 (r - 2M) 2 A" - rM(r - 2M)A' + 2M(2r - 3M)A (130) 
+ (2r - 5M)(r - 2M) 2 C - 6(r - 2M) 2 C + (2r - 3M)(r - 2M) 2 F' 
+6(r - 2M) 2 F + 6r(r - 2M) 3 H" + 6(2r - 3M)(r - 2M) 2 H'. 

We first solve ffT23l) for A, 

A = ^ [C + F + r{r- 2M)H" + 2(r - 2M)iy'] , (131) 

and then solve ^ - (r - 2M) 2 R { }} - R$ for F, 

F = - \r{r - 2M) 2 H"' - 2(r - 2M)(r + 2M)iJ" - 2(5r - 7M)if' - MC'\ . (132) 
6 
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The remaining equations then need to be solved for H' . By combining equations, we 
can obtain a second-order inhomogeneous equation for H', 

- 2r(r + M)(r - 2M) 2 H"' - 2(4r 2 + 3rM - 4M 2 )(r - 2M)H" (133) 
+2(4r 2 - 13rM + AM 2 )H' = M[(r - 2M)C + QC}. 

It's useful to note at this point that if M = 0, we have a solution with F = H' = 
and A = C = —fi 2 r 2 , which is precisely our original plane wave. 
The general solution of (11331) is 

_ /j, 2 r 2 - 2M 2 c [-QrM(r + M) + 4M 3 + (6rM 2 - 3i? 3 ) ln(l - 2M/r)] 

- yCr- )+Ci - _ 2M +c 2 — - — 

(134) 

This then satisfies all of the equations. To get a solution which is both regular and has 
the correct asymptotics, i.e. has A — > — /i 2 r 2 at large r, we need to take ci = — |/i 2 
and C2 = 0. We find 

3 r - 2M ' v ; 



and 

16 , r2 M 3 
— M - 2— 

3 r 



A = -ji 



r 2 - ArM + —M z - 2— 



F _ 2/i 2 M3r 2 -9rM + 5M 2 

~ 3 r - 2M ' ^ ^ 



In terms of the alternative combinations A, F, 

M 21 

r-2M+ 

3r 



A = -fi 2 (r - 2M) 



F = fi 2 M(2r - M). (137) 



Thus, this solution satisfies the regularity conditions at the horizon. Regularity of 
the original functions a(r), f(r),g(r), h(r) at r = 2M further requires us to choose 

2u 2 M 2 

g(2M) + 3h{2M) = — " . (138) 

3 

We now match the near horizon and intermediate region solutions in the interme- 
diate region yU -1 3> r ^> M, where both approximations are valid. The contribution 
from the black string background is 

ds 2 NR>BG ~ -(1 - ^~)dt 2 + (1 + ^~)dr 2 + r 2 {d6 2 + sin 2 Odcj) 2 ) + dz 2 . (139) 

We must now find the unknown functions a(r), b(r), c(r), f(r),g(r), h(r) in this region 
to obtain the contribution from the perturbation. In addition to the solutions 



(11291) and (11361) we must make a choice of gauge. We choose g + 3h = —Mfj, 2 r 
in order to make the rr-component of the perturbation vanish, matching our gauge 
choice in the intermediate region solution. We find, keeping just the terms up O(M) 
and OGu 2 ), 

a(r) w -/i 2 (r 2 - 4Mr), b(r) w -yU 2 (r 2 - 3Mr), c(r) w -/i 2 (r 2 - 2Mr), (140) 
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/(r) « 0, #(r) « -M/i 2 r, h(r) « 0. (141) 
Hence the near region perturbation is 

ds 2 NR,p ~ (a(l - 3cos 2 #) + / 3sin 2 fl(cos 2 0- sin 2 0)) x (142) 
(-fi 2 r 2 (dt + dzf + M/i 2 r(4dt 2 + Qdtdz + 2dz 2 - r 2 (d9 2 + sin 2 0# 2 ))). 

In the intermediate region the plane wave background is, 

ds 2 IRtBG = -dt 2 + dz 2 + dr 2 + r 2 (d6 2 + sm 2 6d(f) 2 ) (143) 
-// 2 r 2 (a(l - 3 cos 2 9) + j3 sin 2 #(cos 2 <p - sin 2 <j>)){dt + dz) 2 . 

From section 12.31 the perturbation due to the black string is 

ds) RP = ^-dt 2 + ^-dr 2 + M/i 2 r(a(l - 3 cos 2 9) + f3 sin 2 6>(cos 2 - sin 2 0)) x 
(Adt 2 + Qdtdz + 2dz 2 - r 2 (d^ 2 + sin 2 9d(f) 2 )). (144) 

Thus the solution constructed in the near region dsf^j^ — dsf^j^ R q -\- dsf^j^ p agrees 
with the solution constructed in the intermediate region ds 2 R = ds 2 R BG + ds 2 IR P to 
the relevant order. This gives us an approximate solution describing a black string 
in a plane wave, valid when the size of the black string is small compared to the 
curvature scale of the wave, r + 

As in [13], the perturbation does not affect the thermodynamic properties of the 
black hole at this order. The area of the horizon cannot be affected at this order 
because the perturbation is entirely in an I = 2 mode, which deforms the shape of the 
S 2 but does not change its area. The temperature cannot be affected because it is 
constant over the horizon. Since the perturbation is an I = 2 mode, it will vanish at 
some point on the horizon, whence the temperature at that point must be unaffected, 
and since it is constant, it must be unchanged over the whole horizon. 



4 Conclusions 

In this paper, we have attempted to construct solutions describing black holes and 
black strings in plane wave backgrounds using the matched asymptotic expansion 
method. We have found that it is not possible to construct a regular black hole solu- 
tion. In the approximation where the wave is thought of as a linearised perturbation 
on the black hole solution, we need a non-zero vector part in the spherical harmonic 
decomposition on the sphere, and it is not possible to make this vector part regular 
on the horizon. It would be interesting to have a deeper physical understanding of 
this failure of regularity. One might think that this is simply saying that the plane 
wave is exerting a force on the black hole, so no stationary solution exists. However, 
we do not believe this is the correct interpretation of our result. The black hole was 
chosen to follow a geodesic in the plane wave background, so there is no force on it 
at leading order. Finite size effects can be analysed in the asymptotic region using 
the classical effective field theory approach of [23 ED] • In this approach, the work 
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done by such finite size terms involves derivatives of the long wavelength background 
fields along the black hole world-line. Since our world-line is chosen to be an orbit of 
the isometries of the background, the work done will vanish. Thus, we would have 
expected the background to simply produce some deformation of the horizon. 

The regularity problem seems to be simply an inconsistency between the symmetry 
structure of the black hole and the plane wave. In four dimensions, the problem is that 
the solution will not be axisymmetric, so there cannot be a regular black hole solution 
as all stationary four-dimensional black holes are required to be axisymmetric |25j . 
In higher dimensions, however, stationary axisymmetric solutions describing black 
holes in plane waves could in principle exist, and the fact that our solutions are never 
regular is somewhat mysterious. Further exploration of this issue is an interesting 
project for the future. 

The importance of this problem is reinforced by the fact that the failure of regular- 
ity here is a counter-example to the assumption in [16] that satisfying the blackfold 
equations implies horizon regularity. Understanding this issue in a more general 
context is clearly important for the blackfolds program [151 US]; i n considering the 
embedding of black branes in arbitrary backgrounds, we need to understand when 
the resulting deformation of the near-horizon region will preserve the regularity of the 
event horizon. Clearly we must require that the embedding of the blackfold in the 
background spacetime preserves enough symmetry to satisfy the rigidity theorems of 
[251 [27]. Ovrr higher- dimensional examples indicate that this is a necessary but not a 
sufficient condition. Identifying sufficient conditions is an important general problem. 

We successfully constructed an approximate solution describing a black string in 
an asymptotically vacuum plane wave background in 5 dimensions. It would clearly 
be interesting to extend this work to find black string solutions in backgrounds which 
asymptote to maximally supersymmetric plane waves. It should be straightforward 
to extend our calculation to this case. 

Our analysis has also led to an interesting general result; the effect of localised 
objects in a plane wave background is not small, even far from the source. The usual 
1/r falloff associated with a localised object in d + 1 spatial dimensions is offset 
by the fi 2 r 2 factors coming from the plane wave background. As a result, we find 
that the "perturbation" due to the source is larger than the background metric at 
sufficiently large r. This leads us to believe that these solutions should not be thought 
of as "asymptotically plane wave" spacetimes. 

A definition of "asymptotically plane wave" was proposed in our earlier work |17j . 
which allows the construction of a well-behaved action principle. This still seems a 
useful definition. However, from the present results it seems that the phase space as- 
sociated with those boundary conditions will not include solutions describing localised 
sources in a plane wave background, so it may not admit many physically interesting 
solutions. Understanding the space of asymptotically plane wave spacetimes is clearly 
important for attempts to construct a direct holographic duality directly for plane 
waves, so we would like to understand this issue better. 

Similar problems have arisen in AdS2 spacetimes [18J, where there are no finite- 
energy asymptotically AdS2 geometries, and in the study of near-horizon extremal 
Kerr solutions (NHEK) [T^l [201 EI]> where the space of metrics which are asymp- 
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totically NHEK consists only of the NHEK solution and solutions obtained from it 
by diffeomorphisms. It is interesting to note that the plane waves, like AdS2, have 
a one-dimensional boundary [3TJ E2]- Perhaps the problem is that there is in some 
sense "not enough space" near infinity to have interesting asymptotically plane wave 
solutions. It would be interesting to carry out a general analysis for asymptotically 
plane wave solutions along the lines of that in [20j 121] ■ We leave this as a project for 
the future. 
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A Appendix 

In this appendix we prove some harmonic identities needed for our analysis of black 
holes in the near horizon region. Definitions are given in section 13.11 We want to 
show that: 



• D % D 3 ViA 







Proof: 



oc D'DWiVj + D'DWjVi 

= [D\ D j ]DiVj + m&DiVj 

= -R k i^D k V 3 - R k / j DiV k - 2k 2 v D^ 

= RUDtVj-lFDiVk 

= 0. 



(145) 



Proof: 



D l D 3 S %1 



k 1 2 

-ja^W D^DjS + ^DWiDWjS + ^DWjS 
-±- 2 D\R k / l S) + f jS 
^D\R\D k S) + 



(146) 



for S 2 , R, 



■ij Tij SO 



D l D 3 Sa 




(147) 
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. D% j = -^(k 2 -2)D j S, 
Proof: 



£>% = -^D i D i DjS + ^DjS (148) 



1 



2k2 (P-2) Dj S 



Proof: 



D% = -^(DW^ + D'D^ (149) 
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